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Mid-Infrared Properties of Powerful Radio Sources 
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ABSTRACT 

We have measured mid-infrared radiation from an orientation-unbiased sample of 3CRR galax- 
ies and quasars at redshifts 0.4 < z < 1.2 with the IRS and MIPS instruments on the Spitzer 
Space Telescope. Powerful emission (L 2 4^ m > 10 22 4 WHz~ 1 sr _1 ) was detected from all but 
one of the sources. We fit the Spitzer data as well as other measurements from the literature 
with synchrotron and dust components. The IRS data provide powerful constraints on the fits. 
At 15 /im, quasars are typically four times brighter than radio galaxies with the same isotropic 
radio power. Based on our fits, half of this difference can be attributed to the presence of non- 
thermal emission in the quasars but not the radio galaxies. The other half is consistent with dust 
absorption in the radio galaxies but not the quasars. Fitted optical depths are anti-correlated 
with core dominance, from which we infer an equatorial distribution of dust around the central 
engine. The median optical depth at 9.7 /im for objects with core-dominance factor R > 10~ 2 is 
« 0.4; for objects with R < 10~ 2 , it is w 1.1. We have thus addressed a long-standing question 
in the unification of FRII quasars and galaxies: quasars are more luminous in the mid-infrared 
than galaxies because of a combination of Doppler-boosted synchrotron emission in quasars and 
extinction in galaxies, both orientation-dependent effects. 

Subject headings: active nuclei radio sources — quasars — radio galaxies 



1. INTRODUCTION 

From an observational perspective, active 
galactic nuclei are diverse. The existence of rela- 
tivistic jets in some AGN, with assoc iated Doppler 
boosting (jBlandford fe Reesl Il978f ). guarantees 
that some of this diversity of appearance must 
be due to orientation. It is equally clear that not 
all differences can be due to orientation, and that 
there are real physical differences between AGN. 
Much observational and theoretical work over the 
last three decades has been devoted to sorting 
out orientational and physical differences. The 
greatest successes have been achieved in the "uni- 
fied model" of Fanaroff-Ri ley Class II galaxies 
and quasars (|Barthe]|[l989l) . and in the unifica- 
tion of Seyfert 1 and Seyfert 2 galaxies (see e.g., 
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Antonucci fc Millerlll985T ). In both cases a dusty 
torus or disk-like structure obscures optical and 
ultraviolet emission from the accretion disk and 
environs along some lines-of-sight, accounting for 
apparent differences between thes e classes of AGN 
(see e.g.. lUrrv fc Padovani 19951 ). 

The FRII/quasar case is particularly impor- 
tant for the following reason. These objects are 
selected to have low frequency (178 MHz) radio 
luminosity greater than about 10 25 WHz _1 sr -1 
rtFanaroff fc R ilcv 1974). Low frequency emission 
is synchrotron emission from giant, optically thin 
clouds that therefore emit isotropically. The low 
frequency emission of dense and possibly Doppler- 
boosted regions of these sources is suppressed 
by self-absorption, which produces a v 2 5 spec- 
trum. As a result, bolometrically insignificant 
but highly anisotropic emission, which may dom- 
inate the source at higher frequencies, is almost 
always insignificant at a few hundred megahertz 
or less. Moreover, high spatial resolution observa- 
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tions of these sources with interferometers allow 
one to measure and subtract the small contribu- 
tion of potentially anisotropic emitting regions at 
low frequency. Of all the identifying signatures of 
nuclear activity, only low frequency radio emis- 
sion originates far from the nucleus. Only low 
frequency radio emission comes from outside the 
regions where anisotropic emission and obscura- 
tion are potentially important. Low frequency ra- 
dio emission, therefore, provides a unique way to 
obtain an orientation-unbiased sample of AGN. 

Only a small fraction of AGN produces the 
prodigious radio power of FRII galaxies and 
quasars. Nevertheless, the fact that FRIIs can 
be studied in an orientation-unbiased sample, se- 
lected on the basis of low frequency radio emis- 
sion, makes them uniquely valuable in separating 
the effects of orientation from physical differences. 
Visible-light and X-ray selected samples do not 
have this property. 

The mid- and far-infrared properties of these 
powerful radio sources are largely unknown. Their 
space density is so low that only a few (e.g., 3C 405 
= Cygnus A) are at low redshifts. The Spitzer 
Space Telescope ( Werner et al. 2004 ) promised a 
major advance in sensitivity over previous tele- 
scopes, and the capability to measure these ob- 
jects. We therefore undertook observations with 
Spitzer of powerful radio sources. The overall 
goal was straightforward: to measure for the first 
time the mid- and far-infrared emission from an 
orientation-unbiased sample of the powerful radio 
sources. 

More specific goals can also be considered. 
For example, the central proposition of FRII 
galaxy/quasar unification is that FRII galaxies 
are just quasars seen from an angle in which 
the optical-to-uv emission of the central re- 
gion is intercepted by the opaque, dusty torus, 
which re-emits this radiation at longer wave- 
lengths. Conversely, quasars are just FRII 
galaxies seen closer to the jet axis. Mod- 
els suggest that the dusty torus becomes op- 
tically t hin at wavelengths in the far-infrared 



(see e .g.. iPier k Krolik 11992 



Schartmann et all 120051: IFritz et alj 120061 ), sug- 
gesting a direct test of FRII galaxy /quasar 
orientation-based unification. If radio-lobe lu- 
minosity is an indicator (albeit time-averaged) 
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of central engine energy, then FRII galaxies and 
quasars with comparable isotropic radio-lobe lu- 
minosity should exhibit "comparable" isotropic 
far-infrared luminosity. Comparable does not 
necessarily mean equal, but there should be no 
FRIIs with low nuclear infrared luminosity at 
wavelengths where the emission is optically thin. 

This test has been performed by other inves- 
tigators using infrared measurements from the 
Infrared Astronomical Satellite (I RAS) and the 
Infra red Space Observa t ory (I SO; iKessler et al 



1996). iHeckman etaD {1992) compared IRAS 



measurements of 3CR quasars and radio galax- 
ies (z > 0.3) selected on the basis of their 178 
MHz radio flux density. They found that the rest- 
frame 6-50 /xm emission is on average four times 
stronger in the quasars than in the radio galaxies, 
suggesting that either the quasars are intrinsi- 
cally stronger sources of mid-far infrared emission 
than quasars, or that this emission is anisotropic. 
The latter requires large optical depths in the 
far-infrared, possibly requiring implausibly large 
dust masses. The former possibility is problematic 
for pure orientation-based unification of quasars 
and radio galaxies, but may be due to a significant 
non-thermal contrib ution in th i s wave length range 
from the radio jet. iHes et al.l (| 19951 ) re-analyzed 
the IRAS 60 /im data for 3CR quasars and ra- 
dio galaxies in the redshift range 0.3 < z < 0.8 . 
They confirm the finding of lHeckman et al" ( 1992h 
that 3CR quasars are systematically brighter than 
radio galaxies in the mid-far infrared and con- 
clude that a beamed 60 (im component may ac- 
count for this d ifferen ce. In subsequent work, 
Hoekstra et al.l (|1997l ) find that the observed 
IRAS 60 /im infrared-to-radio flux-density ratios 
are consistent with an orientation-based model in- 
corporating beamed emission, but conclude that 
other processes such as torus optical depth also 
contribute. 

The ISO photometry for four pairs of 3CR 
quasars and radio galaxies mat ched in redshift and 



radio power was compared by Ivan Bemmel et al 
(2000). Although the non-thermal contribution 



was estimated at less than 2%, a systematic in- 
frared excess was found for the quasars with re- 
spect to the galaxies, suggesting the presence of 
optically thick dust. However the small sample 
size di d now allow any firm con clusions to be 
drawn. iMeisenheimer et"aL ( 2001 ) compared ISO 
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photometry for ten pairs of 3C quasars and radio 
galaxies, and found no difference in the distribu- 
tion of their (radio- lobe normalize d) far-infrared 
luminosities. lAndreani et al.l (|2002f ) obtained mil- 
limeter data and ISO photometry for a sample of 
3C radio galaxies and quasars (selected on the ba- 
sis of visibility to ISO) and found that the quasar 
composite spectrum is elevated by over a factor 
of three compared with the radio galaxy compos- 
ite spectrum in the millimeter region, which the 
authors attribute to a beamed non- thermal com- 
ponent. ISiebenmorgen et al. ( 2004f ) apply radia- 
tive transfer models to all 68 3CR sources imaged 
by ISO. The model spectra of Type 1 objects (i.e. 
viewed at inclinations close to the radio axis) peak 
around 40 /im while those of Type 2 (i.e. viewed 
at inclinations away from the radio axis) peak at 
wavelengths over 100 /Ltm. This is interpreted as 
being in su p port of orientation-based unification. 



Haas et al.l (|2004J ) examined the spectral energy 
distributions of 3CR quasars and radio galaxies 
based on ISO photometry, millimeter observations 
using MAMBO and supplementary data from the 
Nasa/IPAC Extragalactic Database (NED). The 
sample was selected on the basis of good visibility 
to ISO. At rest-frame 70 /xm, the authors found 
that the (radio-lobe normalized) thermal dust lu- 
minosity of quasars and galaxies had similar distri- 
butions, in line with orientation-based unification. 

It is clear from previous work that beamed syn- 
chrotron emission and dust extinction modulate 
the mid-infrared emission of quasars and radio 
galaxies to some degree. However, it has not been 
possible to assess quantitatively whether a combi- 
nation of these two effects is sufficient to unify an 
orientation-unbiased sample of quasars and radio 
galaxies. The increased photometric sensitivity 
and spectroscopic data provided by Spitzer pro- 
vide additional constraints in the infrared, allow- 
ing us to perform a spectral decomposition of the 
thermal and non-thermal components of emission 
for the most powerful radio sources and also to 
estimate the silicate optical depth. 

In this paper, we describe the selection method- 
ology (§HJ), Spitzer observations (§[3]) and data re- 
duction (§HJ for our sample of powerful 3C radio 
sources. The photometric and spectroscopic re- 
sults are presented (§ 15. ip and the spectral energy 
distributions (SED) are fitted with a combination 
of synchrotron and dust components (§ 15. 2| . We 



then explore what these results can tell us about 
the physical conditions (§ 16. ip and the origin of 
the infrared emission (§ 16. 2[ § 16. 31) in these ob- 
jects. Finally, the comparative infrared luminosity 
of quasars and galaxies is discussed (§ 16. 4p . 

2. THE SAMPLE 

We require a sample selected at low frequency, 
with 1/178 MHz > 10 26 WHz -1 sr -1 and with a rea- 
sonable balance between FRII radio galaxies and 
quasar^]. 

The sample of Barthell (1989) provides the 
ideal starting point. It contains the 50 sources 
in the complet e low frequency 3CRR catalogue 



of lLaing et all (|1983I ) with 0.5 < z < 1.0. 
All have emitted radio luminosity Ln$ MHz > 
10 26 WHz-'sr 1 . Of the 50, 33 are identified 
with galaxies and 17 with quasars, with classifi- 
cations based on the optical type designated by 



Laing et al.l ( 19831 ). On the basis of radio mor- 
phology, all sources in the sample are FRII, as 
expected from the high radio luminosity. 

To reduce the observing time required, we re- 
duced the sample to 33 objects based on ecliptic 
latitude. Sources at high ecliptic latitude are eas- 
ier to schedule with Spitzer. We also added one 
source to the sample, 3C200, at z = 0.458, be- 
cause we already had a 16 ks Chandra observation 
of it. We emphasize that neither of these changes 
introduces a bias in orientation in the sample. 

The Spitzer sample (hereafter, the "primary" 
sample) then consists of 16 quasars and 18 galax- 
ies. The general properties of the sample arc sum- 
marised in Table [TJ Figure [T] shows the radio 
lobe luminosity of the sample at an emitted fre- 
quency of 178 MHz as a function of redshift. The 
galaxies and quasars span similar ranges in red- 
shift and radio-lobe luminosity. The mean (me- 
dian) redshift for quasars in the primary sam- 
ple is 0.74 (0.73) and for galaxies is 0.79 (0.78). 
The mean (median) radio luminosity, £i78MHz> 
is 5.0(3.8) x 10 27 WHz _1 sr _1 for quasars and is 
3.2(2.6) x 10 27 WHz- 1 sr" 1 for galaxies. The fit to 



1 Since the low frequency radio emission is isotropic, one 
would be justified in multiplying by An and writing the full 
luminosity per hertz; but since luminosities at other fre- 
quencies are not necessarily isotropic, we will generally use 
units of WHz _1 sr _1 . We express even the isotropic low 
frequency radio value in the same way for easy comparison. 
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Fig. 1. — Luminosity (WHz -1 sr -1 ) at rest-frame 
178 MHz vs. redshift. Objects from our sample are 
labeled as follows; closed circles: quasars — closed 
triangles: steep-spectrum compact (SSC) quasars 
— closed squares: superluminal quasars — closed 
diamonds: SSC superluminal quasars — open cir- 
cles: galaxies — open triangles: SSC galaxies. 

the synchrotron SED used to interpolate the ob- 
served low-frequency flux density to the emitted 
frame is described in § 15.21 The sample includes 
five steep-spectrum compact (SSC) quasars and 
one SSC galaxy. These objects are characterised 
by small projected linear sizes (< 15 kpc) and 
convex radio spectra peaking around 100 MHz, 
steepening t o a < —0.5 (S v oc v a ) at higher 



frequencies ( O'Dea 19981 ). The sample also in- 



cludes seven quasars with reported superluminal 
motion, 3C138, 3C 147, 3C216, 3C263, 3C334, 
3C336 and 3C380. 

For comparison, we will show results from ISO 
observati ons of six 3CRR o bjects from the compi- 
lation of Haa s et al. (2004). Sp itzer ob servations 



of thirteen 3CRR objects from IShi et all (|2005l ). 



Spitzer observations of ten 3CRR radio galaxies 
from Birkinshaw et al. (2006, in preparation) and 
Spitz er observat i ons o f nineteen 3CRR galaxies 
from I Ogle et al.l (2006). Table [2] summarises the 
properties of these 48 additional objects which 
comprise the supplementary sample. Together, 
the primary and supplementary samples comprise 
82 3CRR objects. 



3. OBSERVATIONS 

All sources in the primary sample were sched- 
uled for observation with the Long-Low mod- 
ule of the Spitze r Infrared Spectrograph (IRS; 
Houck et al ■112004)) resulting in spectra in the range 
15-37 /im, as well a s with the Multib and Imaging 
Photometer (MIPS; iRieke et al1l2004l ) in photom- 
etry mode at 24, 70, and 160 /im. Of the 34 obser- 
vations with IRS, four (3C200, 3C 220.3, 3C272 
and 3C 289) peaked-up incorrectly, either on tran- 
sient pixels or on a brighter object in the field. 
The galaxy 3C 427.1 was not detected with IRS 
and 3C6.1 was detected in LL1 only. All objects 
in the sample are unresolved by the Spitzer instru- 
ments. 

4. DATA REDUCTION 

The IRS spectra were extracted from the 
basic calibrated data (BCD) images provided 
by the Spitzer Science Center, using the Spec- 
troscopy Modeling and Analysis Tool (Smart; 
Higdon et al.l2 004) . The images were background- 
subtracted by forming difference images from the 
BCD data at the two nodding positions. Spec- 
tra were then extracted using narrow extraction 
apertures to maximise the signal-to-noise ratio, 
and calibrated using a standard IRS Long-Low 
stellar calibrator, HD 173511, correcting the flux 
calibration for the width of the narrow extrac- 
tion aperture. The spectra were rescaled to the 
MIPS 24 /im point in order to provide a consistent 
flux calibration between MIPS and IRS measure- 
ments. The MIPS and IRS observations were not 
contemporaneous; if the source is variable, this 
cross-calibration introduces an error. 

Mosaics of the MIPS 24 and 70 fim BCD images 
were produced using the Mos aicking and Point 
Sourc e Extraction tool, Mopex (jMakovoz fc Khan 
2005), while for 160 fim, the post-BCD mosaics 
provided by the SSC were used. The radius of the 
photometric aperture was set to 15", 15", and 48" 
at 24, 70, and 160 /mi, respectively. Variations in 
the background on the scale of the aperture were 
used to estimate the uncertainty in the extracted 
flux densities. 
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Table 1 

Summary of primary 3CRR sample properties. 



Source 


ID" 




z 


l°g £-178 MHz 

[WHz -1 si- -1 ] 


fiapp 


c 


3C 138 


Q(SSC)(SL) 





759 


27.43 


2.1h~ 




3C 147 


Q(SSC)(SL) 





545 


27.48 


1.3/1" 


1 


3C 175 


Q 





768 


27.60 






3C 196 


Q 





871 


28.25 






3C207 


Q 





684 


27.36 






3C216 


Q(SL) 





67 


27.48 


2.6h~ 


l 


3C 254 


Q 





734 


27.61 






3C263 


Q(SL) 





646 


27.32 


1.3/1" 


l 


3C 275.1 


Q 





557 


27.14 






3C 286 


Q 





849 


27.67 






3C 309.1 


Q(SSC) 





904 


27.76 






3C 334 


Q(SL) 





555 


27.03 


2.0h~ 


l 


3C 336 


Q(SL) 





927 


27.64 


2.6h~ 


l 


3C 343 


Q(SSC) 





988 


27.40 






3C 380 


Q(SSC)(SL) 





691 


27.99 


8.0h~ 


l 


3C6.1 


G 





84 


27.54 






3C 22 


G 





937 


27.63 






3C 184 


G 





994 


27.65 






3C200 


G 





458 


27.86 






3C 220.1 


G 





62 


27.26 






3C 220.3 


G 





685 


27.39 






3C 268.1 


G 





9737 


27.80 






3C 272 


G 





944 


27.46 






3C280 


G 





9975 


28.94 






3C289 


G 





9674 


27.65 






3C 292 


G 





713 


27.16 






3C 325 


G 


1 


135 d 


27.88 






3C 330 


G 





55 


27.41 






3C 337 


G 





635 


27.19 






3C 340 


G 





7754 


27.30 






3C 343.1 


G(SSC) 





75 


26.13 






3C 352 


G 





8057 


27.50 






3C 427.1 


G 





572 


27.48 






3C 441 


G 





707 


27.35 







a G — galaxy, Q — quasar, SSC — steep spectrum compact and 
SL — supcrluminal source. 

b scc §ora 

c Data on superluminal motion taken from [Zensus et alj ([1987^ , 
IVenturi et al.l U993|). | Mantovani et al.l I1998B . [Shen et al.l 0001), 
and lHough et al.l I2002B . 

d We use t h e upd ated redshift of z — 1.135 for 3C 325 from 
IGrimes et alj ||2005I 1. 
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Table 2 

Summary of supplementary 3CRR sample properties. 









log i e (178 MHz) 




Source 


ID 




[\VHz~ 1 sr -1 ] 


Reference 





Q 












3C 48 


'■J 


|| 


■<"-!' 


%~ 


29 




3C 109 


Q 





.3056 


26 


.73 


1 


3C 249.1 


Q 





.3110 


26 


44 


2 


3C 351 


Q 





.3710 


26 


.71 


2 


3C 454.3 


Q(SL) 





.8590 


27 


.35 


1 


3C 28 


G 





.1952 


26 


.19 


4 


3C 33 


G 





.0595 


25 


.60 


4 


3C 33.1 


G 





.1810 


25 


.99 


1 


3C 55 


G 





.7350 


27 


.67 


4 


3C61.1 


G 





.1860 


26 


.41 


2 


3C66B 


G 





.0215 


24 


34 


3 


3C 83. IB 


G 





.0255 


2-1 


.53 


3 


3C 123 


G 





.2177 


27 


.33 


4 


3C 153 


G 





.2769 


26 


47 


4 


3C 171 


G 





.2384 


26 


.45 


2 


3C 172 


G 





.5191 


27 


.11 


4 


3C 173.1 


G 





.2920 


26 


55 


2 


3C 192 


G 





.0598 


25 


.19 


3 


3C 219 


G 





.1744 


26 


.47 


2 


3C 234 


G 





.1848 


26 


.41 


1 


3C 236 


G 





.0989 


25 


47 


2 


3C 244.1 


G 





.4280 


27 


.04 


4 


3C 263.1 


G 





.8240 


27 


.67 


4 


3C 265 


G 





.8108 


27 


.71 


4 


3C 274.1 


G 





.4220 


26 


94 


4 


3C284 


G 





.2394 


26 


22 


2 


3C 285 


G 





.0794 


25 


.18 


3 


3C 288 


G 





.2460 


26 


.46 


4 


3C 300 


G 





.2720 


26 


.53 


4 


3C 310 


G 





.0540 


25 


.52 


4 


3C 315 


G 





.1083 


25 


.65 


2 


3C 319 


G 





.1920 


26 


.14 


4 


3C 321 


G 





.0960 


25 


.42 


3 


3C 326 


G 





.0895 


25 


.54 


4 


3C 349 


G 





.2050 


26 


.13 


2 


3C 381 


G 





.1605 


25 


.99 


2 


3C 382 


G 





.0578 


25 


.13 


3 


3C 386 


G 





.0177 


24 


.16 


3 


3C 388 


G 





.0908 


25 


.63 


2 


3C 390.3 


G 





.0569 


25 


.50 


1 


3C401 


G 





.2010 


26 


.30 


4 


3C433 


G 





.1016 


26 


.10 


4 


3C436 


G 





.2145 


26 


.31 


4 


3C438 


G 





.2900 


27 


00 


4 


3C442A 


G 





.0270 


24 


.37 


3 


3C452 


G 





.0811 


25 


88 


3 


3C465 


G 





.0293 


24 


.81 


3 



es. — (1)[H 
t et al. (2006 



aas et al.l 120041 1. ( 21 IHhi et al.l 12005ft . 
in preparation), (4) |0gle ef^TTIiioM) 



G 



5. RESULTS 

This section presents the infrared measure- 
ments for the primary and supplementary sam- 
ples. The infrared luminositie^ of the sample are 
characterised at rest-frame wavelengths of 15 and 
30 /im. 

A rest wavelength of 15 /im was chosen because 
it lies between the bro ad silicate features around 
9.7 and 18 pm (s ee e.g., Siebenmorgen et al. 2005; 
Hao et all [20051) and avoids the [Ne v]A14.3/im 
and [Ne III]A15.5 /im lines. All objects in the 
primary sample are at redshifts such that rest- 
frame 15 /im lies within IRS Long-Low. For four 
objects without IRS spectra (3C 200, 3C 220.3, 
3C272 and 3C289), the rest-frame 15 /im flux 
density was estimated by assuming a power-law 
SED between the MIPS 24 and 70 /im measure- 
ments. The galaxy 3C 200 was detected in the 
MIPS 24 /im channel but not at 70 /im, so the de- 
rived 15 and 30 /im flux densities are considered 
lower and upper limits respectively. The galaxy 
3C 427.1 was not detected using MIPS or IRS. 
However, since rest-frame 15 /im corresponds to 
an observed wavelength of ~ 24 /im for this object, 
we use the observed MIPS 24 /im measurement to 
place an upper limit on the rest-frame 15 /im flux 
density. This galaxy is of interest since it is the 
faintest 15 /im object in the primary sample. The 
rest-frame 30 /im flux density was determined by 
assuming a power-law SED to interpolate between 
the observed MIPS 24 and 70 /im flux densities. 
Where the 70 /im measurement resulted in an up- 
per limit, the interpolated 30 /im flux density is 
also considered an upper limit. The 178 MHz flux 
densities are corrected to the rest-frame using a 
fit to the radio synchrotron spectra as described 
in § 15.21 Table [7] in Appendix A gives the charac- 
teristic 15 and 30 /im flux density and luminosity 
for objects in the primary sample as well as the 
rest-frame 178 MHz luminosity. 

For the supplementary sample, power-law SEDs 
are also assumed in order to interpolate the avail- 
able photometry (e.g., ISO, MIPS, or IRAC) to 
rest-frame wavelengths of 15 and 30 /im. Only 
15 /im data were available for the object s in th e 



supplementary sample from Ogle et al. (2006). 



Extrapolation to rest- frame 178 MHz was achieved 



2 We used H = 70 km s" 1 Mpc -1 , fl M = 0.3 and f2 A = 0.7. 



using the observed 178 MHz flux density and the 
178-750 MHz spectral index. 

5.1. MIPS and IRS measurements 

Table 3 gives the MIPS 24, 70 and 160 /im flux 
density measurements for the primary sample. All 
objects in the primary sample were detected at 
24 /im except for the galaxy 3C 427.1. At 70 /im, 
6 out of 19 galaxies and all the quasars were de- 
tected. At 160 /im, there was just a single detec- 
tion, the galaxy 3C 220.3. 

Figure [2] shows the IRS Long-Low spectra for 
the quasars (left panel) and galaxies (right panel) 
plotted in (arbitrarily scaled) units of luminosity, 
vL v , at the emitted wavelength, arranged in order 
of decreasing total dust luminosity (see § 15. 2j) from 
top to bottom. This figure allows the general mor- 
phology and features to be compared (the spectra 
are presented in greater detail with error estimates 
in Figure [L3|) . Table 4 gives the emission line flux 
densities for the objects in the primary sample. 

The [Ne iii]A15.5 /im line is present at some 
level in most of the galaxies and [Ne v]A14.3/im 
is detected in two galaxies, 3C 184 and 3C330. 
The [Ne iii]A15.5 /im line can also be seen in 
many of the quasars in the sample, as well as 
[Neii]A12.8/im, [Nev]A14.3/im, [S iii]A18.7/iin, 
and [Siv]A10.5/im. [NeVl]A7.6/iin is detected in 
two quasars, 3C 309.1 and 3C343. 

The galaxy continua are generally fainter in 
the mid-infrared than those of the quasars and 
so the signal-to-noise ratio is typically lower for 
the galaxies. The presence of silicates in emission 
is evident in some quasars (e.g., 3C263, 3C 175, 
3C380). Silicate absorption is not obvious from 
a visual inspection of the spectra, except for the 
quasar 3C 343 and the galaxy 3C 184. 

5.2. Spectral Fitting 

In this section we fit the IRS and MIPS data, as 
well as additional photometry from the NED, with 
models for the thermal and non-thermal emission, 
in order to separate their respective contributions 
to the mid-infrared. Figure [3] shows a cartoon 
of the physical origin of the broadband spectral 
components used in the fits: (i) synchrotron emis- 
sion from the radio lobes, (ii) synchrotron emission 
from the radio jet, and (iii) thermal emission from 
dust. Four combinations of components are fitted 
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Fig. 2. — This plot shows the spectra for all the quasars (left) and galaxies (right) in the sample, so that 
the general morphology and spectral features may be compared. The spectra are arranged in order of 
decreasing total dust luminosity (see Table 7). The data for the galaxies have been smoothed over four 
spectral resolution elements. LL2 was not detected for the galaxy 3C6.1. For clarity, the errors on the 
spectra are not shown (see Figure fT3l instead) . 
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Table 3 

mips photometry for the primary sample. 







C24^m 










Name 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


3C6.1 


2.8 


0.3 


<13.8 


4.6 


<45.6 


15.2 


3C 22 


13.4 


1.3 


22.8 


5.6 


<45.6 


15.2 


3C 138 


15.9 


1.6 


24.3 


6.5 


<45.6 


15.2 


3C 147 


22.4 


2.2 


57.2 


12 4 


<88 8 


29 6 


3C 175 


15.5 


1.6 


25.1 


6.4 


<60.0 


20.0 


3C 184 


2.7 


0.3 


10.7 


2.9 


<45.6 


15.2 


3C 196 


18.1 


1.8 


12.2 


3.3 


<95.4 


31.8 


3C 200 


1.3 


0.1 


<10.8 


3.6 


<45.6 


15.2 


3C 207 


11.4 


1.1 


21.5 


6.9 


<77.1 


25.7 


3C216 


20.9 


2.1 


69.5 


14.6 


<45.6 


15.2 


3C 220.1 


2.2 


0.2 


<7.8 


2.6 


<83.4 


27.8 


3C 220.3 


4.8 


0.5 


20.7 


4.8 


136.3 


31.1 


3C 254 


12.4 


1.2 


11.7 


2.5 


<84.9 


28.3 


3C263 


27.0 


2.7 


31.1 


6.5 


<46.5 


15.5 


3C 268.1 


2.4 


0.2 


<6.9 


2.3 


<50.7 


16.9 


3C 272 


2.1 


0.2 


9.8 


2.2 


<85.8 


28.6 


3C 275.1 


8.1 


0.8 


23.6 


7.8 


<45.6 


15.2 


3C 280 


9.3 


0.9 


18.2 


4.0 


<61.2 


20.4 


3C 286 


7.7 


0.8 


29.3 


6.6 


<45.6 


15.2 


3C 289 


3.7 


0.4 


17.4 


4.0 


<81.3 


27.1 


3C 292 


1.9 


0.2 


<8.7 


2.9 


<83.4 


27.8 


3C 309.1 


17.8 


1.8 


36.1 


8.0 


<45.6 


15.2 


3C 325 


3.4 


0.3 


<7.2 


2.4 


<45.6 


15.2 


3C 330 


7.5 


0.7 


<13.5 


4.5 


<45.6 


15.2 


3C 334 


35.0 


3.5 


65.3 


13.6 


<45.6 


15.2 


3C 336 


4.7 


0.5 


12.0 


3.8 


<46.5 


15.5 


3C 337 


1.0 


0.1 


<7.8 


2.6 


<47.7 


15.9 


3C 340 


2.4 


0.2 


<9.3 


3.1 


<68.1 


22.7 


3C 343 


7.8 


0.8 


41.2 


10.8 


<70.8 


23.6 


3C 343.1 


1.3 


0.1 


<9.0 


3.0 


<86.4 


28.8 


3C 352 


1.2 


0.1 


<6.3 


2.1 


<47.4 


15.8 


3C 380 


28.1 


2.8 


34.6 


8.4 


<60.0 


20.0 


3C 427.1 


<0.4 


0.1 


<7.8 


2.6 


<45.6 


15.2 


3C441 


2.4 


0.2 


<4.8 


1.6 


<60.9 


20.3 



Note. — The uncertainties include an absolute flux calibration uncertainty 
of 10% at 24 fim and 20% at 70 and 160 fim. The upper limits are 3<r; 
detections arc in bold font. 
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Table 4 

Emission line flux densities for the primary sample. 





[No n] 




[No ill] 


[No v] 




[Ne vi] 




[Sin] 


[Siv] 




12.8 (im 


15.5 |im 


14.3 fJ,m 




7.6 |im 


18.7/im 


10.5 |im 


Name 


21.6 eV 


41.0cV 


97.1oV 


126. 2oV 


23.3 eV 


34.8 eV 


3C6.1 


< 8.6 




6.3 ± 1.3 


< 9.2 






< 


15.5 




3C 22 


< 4.2 




3.4 ± 0.8 


< 5.1 


< 


5.9 






5.4 ± 1.5 


3C138 


< 9.7 


< 


5.5 


< 8.0 








5.3 ± 0.8 


< 8.8 


3C147 


< 10.5 




8.3 ± 0.7 


< 6.3 








5.3 ± 0.5 


6.6 ± 0.3 


3C184 


1.9 ± 0.4 




6.5 ± 1.1 


1.7 ± 0.5 


< 


4.4 






< 5.7 


3C196 


< 6.4 




4.6 ± 0.4 


3.0 ± 0.4 










< 9.8 


3C216 


< 6.7 


< 


3.8 


< 6.1 






< 


5.2 


< 3.8 


3C 220.1 


< 3.8 


< 


2.1 


< 2.1 






< 


2.1 


< 4.1 


3C 254 


< 8.2 




5.7 ± 1.3 


< 11.9 






< 


18.0 


< 10.3 




< 17.9 




4.7 ± 1.0 


< 9 9 








3.2 ± 0.3 


5 ± 9 


3C 268.1 


< 0.5 


< 


0.9 


< 0.5 


< 


0.6 






< 0.9 


3C 275.1 


< 5.1 




6.7 ± 0.7 


< 5.4 






< 


3.1 


< 5.6 


3C 280 


< 3.6 




4.2 ± 0.5 


< 6.7 


< 


5.2 






< 13.7 


3C286 


< 9.7 




6.9 ± 0.6 


7.5 ± 0.4 






< 


10.4 


7.4 ± 0.9 


3C 292 


< 7.2 


< 


4.6 


< 7.7 






< 


6.1 


< 10.4 


3C 309.1 


< 5.5 




4.3 ± 0.6 


3.3 ± 0.4 




7.2 ± 1.5 






< 10.6 


3C 325 


< 3.7 


< 


4.3 


< 3.3 


< 


4.6 






< 9.7 


3C330 


< 6.1 




4.0 ± 0.2 


4.2 ± 0.6 






< 


2.2 


3.5 ± 0.7 


3C 334 


6.6 ± 0.4 




16.2 ± 1.3 


11.8 ± 1.7 








3.1 ± 0.6 


11.8 ± 0.6 


3C336 


< 4.9 




1.8 ± 0.1 


< 5.4 


< 


5.0 






< 8.9 


3C340 


< 10.9 


< 


3.2 


< 5.6 






< 


6.3 


< 9.5 


3C 337 


< 3.2 


< 


3.1 


< 3.4 






< 


2.4 




3C 343 


7.0 ± 0.9 




13.9 ± 1.3 


8.3 ± 1.1 




6.2 ± 0.9 






7.9 ± 0.6 


3C 343.1 


3.9 ± 0.5 




2.2 ± 0.3 


< 3.8 






< 


5.8 


< 5.2 


3C352 


2.9 ± 0.3 




1.4 ± 0.4 


< 2.7 






< 


4.0 


< 6.4 


3C380 


< 17.1 




4.4 ± 1.1 


< 5.3 






< 


4.3 


4.3 ± 0.6 


3C441 


< 4.9 


< 


5.4 


< 6.2 






< 


5.3 


< 9.7 



Note. — The emission line flux densities are given in units of 10~ 22 W cm . The upper limits are 3cr; detections 
are in bold font. The ionisation potential shown above each column is that needed to create the species from the next 
lower ionisation stage. 
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Fig. 3. — This cartoon indicates the physical origin of the broadband spectral components used to fit the 
spectral energy distributions, using the SED of the quasar 3C 138 as an example (see Fig. H3|). The drawing 
is not to scale — the dust shell has been drawn ps 10 4 times its correct size relative to the extent of the jet. 
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to all objects, and the combination with the best 
reduced x 2 is selected as the best fit. We now de- 
scribe the specific functional forms used for each 
component. 

In order to characterise the synchrotron emis- 
sion from the radio lobes, F}° he , the following 
parabolic function (s ee e.g., Polletta et al.1 liooot 
Andreani et al.ll2002l ) is fit to the data, 



logF' ob » oc -/3(lo gi /-log^ t ) 2 +log[e-^° sc )], 

(1) 

where (3 is a parameter representing the bending 
of the parabola, v t is the frequency at which the 
optical depth of the synchrotron-emitting plasma 
reaches unity and, v l ° hc is the frequency corre- 
sponding to the cut-off energy of the lobe plasma 
energy distribution. The free parameters and con- 
straints for this component are (3, v t > 1 MHz and, 
30 < i/' obc < 3000 GHz. 

A broken power-law is used to represent syn- 
chrotron emission from the jet, i 7 ^ *, viewed 
through a dust s creen, as follows (see e.g., 
Polletta et al]|2000h . 



-{v /v)° 



(2) 

where, r is the optical depth of the dust compo- 
nent (discussed below), a\ and a.2 are the opti- 
cally thick and optically thin spectral indices re- 
spectively, i>q is the frequency at which the broken 
power-law peaks and, u l ° he is the frequency cor- 
responding to the cut-off energy of the jet plasma 
energy distribution. In the fits, the optically thick 
spectral index is fixed at ol\ = 2.5 as expected 
for a homogenous source, and the free parame- 
ters and constraints are r, 3 < uq < 300 GHz, 
6 x 10 13 < t^ ct < 60 x 10 13 Hz, and, a 2 > -1.2. 

The warm dust component is modeled as an 
optically thin spherical shell of graphite and sili- 
cate grains (properties from iDraine &: Led Il984f ) 
surrounding a source with an AGN accretion 
disk spectral energy d istribution (taken from 



Schartmann et alj 120051 ) . Equilibrium tempera- 



tures are calculated for all grains as a function of 
their size and composition. The temperature of 
a dust component is not a well defined quantity 
since the component consists of a distribution of 
grains at different equilibrium temperatures. We 
define the "temperature" of the distribution to 
be the composition-averaged temperature of the 



luminosity-dominating grain size. Grains are dis- 
tributed in si z e with the power-law distribution of 
iMathis et all (|l977l) . The upper grain-size cutoff 
is 0.25 /im and the lower size limit is determined by 
the minimum grain size that has not sublimated 
at the temperature of the distribution. 

For the purposes of producing a simple method 
to characterize spectra, we make the assumption 
that all dust emission in the IRS spectral range 
originates from grains that are directly heated by 
UV photons from the illuminating source. Since 
the extinction in the ultraviolet is more than 100 
times the extinction in the mid-infrared, the dis- 
tance over which all UV photons are absorbed by 
the dust is much less than an optical depth of unity 
in the mid-infrared. We may therefore assume 
that all infrared emission within the spherical shell 
is optically thin. We also allow for the possibility 
that the dust component is partially obscured be- 
hind a screen of cooler dust (defined as dust not 
emitting significantly in the mid- infrared) . There- 
fore, the total model for the emergent dust emis- 
sion is 

^ dust (A) oc£„(A) e - TSC "( A \ (3) 

where E u is the emissivity of the dust distribution 
and r scr is the screen optical depth. 

As an alternative to the screen extinction 
model, we also use a "mixed" extinction scenario, 
where warm dust emitting in the mid-infrared is 
mixed with cooler dust from the shaded regions of 
clouds. This model may be more appropriate for 
a "clumpy" torus. In this case, the dust emission 
model is 



F„ dust (A) oc£„(A) 



lix (A) 



(4) 



where r mlx is the "mixed" optical depth. 

The screen and mixed models represent lim- 
iting cases where the warm and cool dust is ei- 
ther completely separate or completely mixed, and 
should therefore capture the range of possible be- 
haviours. We perform the fits for each object with 
both dust models to find the screen and mixed 
optical depths. 

Our model continuum includes a silicate emis- 
sion feature at 9.7 /im. The optical depths that 
we determine are therefore greater then would be 
determined if we had simply used an interpolated 
smooth continuum. In principle, the extinction 
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values we derive from the fits give a better measure 
of the true extinction to the source, since they in- 
corporate this more realistic continuum. However, 
a full radiative transfer model may provide yet an- 
other optical depth value, since the obscuring dust 
will re-radiate the absorbed power, partially filling 
in the silicate absorption feature. 

Given the poor MIPS detection statistics at 
70 /im and particularly 160 /im, we typically lack 
photometric constraints in the far-infrared. How- 
ever, in order to fit to the far-infrared photome- 
try where available, we also included in the fits a 
cooler dust component, with the temperature al- 
lowed to vary in the range T = 30-80 K. For those 
objects without IRS spectra (3C 200, 3C 220.3, 
3C272, 3C289, and 3C 427.1), we fitted two dust 
components at fixed temperatures T = 80 and 
185 K to the available infrared photometric data. 
(For the galaxy 3C 220.3, which had the only 
160 /im detection, we allowed the temperature of 
the cooler component to be a free parameter). 

For all of the objects in the primary sample with 
IRS spectra, we performed fits to the SED using 
the following combinations of components: 

• warm dust + lobe synchrotron; 

• warm dust + lobe synchrotron + jet syn- 
chrotron; 

• warm dust + lobe synchrotron + cool dust; 

• warm dust + cool dust + lobe synchrotron 
+ jet synchrotron, 

with the free parameters varying simultaneously. 
In general, the combination of components which 
resulted in the best reduced chi-squared was se- 
lected as the best fit. However, an {n + 1)- 
component model is chosen in preference to an 
n-component model only if the inclusion of the 
additional component results in a confidence value 
for the F-test statistic of 99% or better. Figure [T3l 
shows the best fits for the quasars and galaxies in 
the primary sample. Figure fT4l shows how the IRS 
spectra help to constrain the jet synchrotron con- 
tribution by comparing the fits with a jet compo- 
nent replacing cool dust in three quasars. Notes on 
the fits for some individual objects are presented 
in Appendix B. Table[5]gives the parameters of the 
fits as well as the estimated dust luminosity, Ldust , 
and the ratios of thermal to total luminosity at 15 
and 30 /im. The latter constitute the corrections 



which are applied to the characteristic luminosities 
in order to remove the non-thermal contribution. 
The estimated dust luminosity, Ldust, is the inte- 
grated luminosity of the model dust components 
and so does not include any contribution from the 
jet or lobe components. Since we do not model hot 
dust in this work and since we lack constraints in 
the far-infrared, this estimate should be consid- 
ered a lower limit. 

Our spectral fitting methodology is subject to 
the usual problems of fitting simple models of 
complex phenomena to sparse data. The fol- 
lowing caveats are worth pointing out explicitly. 
The spectral energy distributions of the objects 
in our sample are sparsely sampled and we lack 
constraints in the sub-millimeter region in many 
cases. The data used in the fits were not ob- 
tained simultaneously and variations in core flux 
density can reach w 10%. We assume that a sin- 
gle power-law adequately describes the optically 
thin jet synchrotron emission; however, the actual 
synchrotron spectrum is likely to be a complex 
superposition of emission from discrete compact 
regions. We generally lack constraints in the far- 
infrared, so we cannot constrain the existence of 
colder dust. Moreover, the errors on the estimated 
optical depths and non-thermal contributions are 
likely to be dominated by model-dependent uncer- 
tainties that are difficult to quantify. 

Nevertheless, the methodology described above 
introduces no biases into the fitting of compo- 
nents for quasars with respect to galaxies and, 
as discussed § 16.31 below, it is only the quasars in 
the sample which are estimated to have a large 
non-thermal contribution. As a consistency check, 
we also apply an orientation-base d model for the 



beam ed component of emission (jHoekstra et al 



19971 ) to our data in § 16.31 and find that this also 



predicts significant non-thermal contributions for 
the quasars in our sample, suggesting that we 
are not systematically introducing jet components 
where none would be expected. Moreover, all of 
the model components that we use have been well 
established in detailed studies of individual ob- 
jects, and our fitting models are chosen conser- 
vatively. 
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Landscape table to go here 
Table 5: 




10"" 10"" 10" 

[Ne V]X14.3 [W/cm 2 ] 

Fig. 4. — Observed flux density for [Ne v]A14.3 /im 
vs. [Ne ill] A15.5 /im for the primary sample as well 
as other 3CRR sources from lHaas et"al] (|2005l) and 
lOgle et al.l (J2006) (plotted using smaller symbols). 
Only objects with detections in both lines are plot- 
ted. The dashed line (of slope 0.9) is a fit to the 
data for the quasars only. Symbols are as in Fig- 
ure m 



6. DISCUSSION 

We detect powerful emission in the mid-infrared 
for our sample of radio galaxies and quasars. All 
objects in our sample except one were detected 
at 24 //m using MIPS, with luminosities ^ 
10 22 ' 4 WHz -1 sr _1 . In the following section we 
discuss the AGN physical conditions derived from 
the emission line measurements, the origin of the 
measured infrared emission, and the comparative 
luminosity of quasars and radio galaxies. 

6.1. AGN Physical Conditions 

The high-excitation (97.1 eV) [Ne v]A14.3 /im 
line was detected in five quasars (see Table 4); 
it was also detected in two radio galaxies, 3C 330 
and 3C 184, directly confirming the presence of a 
highly luminous ionizing radiation field in these 
sources. We can use the neon emission line ratios 
to probe the physical conditions in the photoion- 
ized region. The [Nev]A14.3/im/[Nelll]A15.5/tm 
ratio responds to the ionization parameter and 
the [Neii]A12.8/ini/[Neill]A15.5/mi ratio is sen- 



c 
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0.5 



log (S /S ) 
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Fig. 5. — This figure shows the distribution of 
the thermal 15 to 30 /im flux density ratio for the 
primary and supplementary samples together (un- 
filled histogram) as well the primary sample alone 
(cross-hatched histogram). The vertical dashed 
lines indicate the approximate division (« 0.17- 
0.21) between AG N and starbu r sts fo und for the 
15/30 /xm ratio bv lBrandl et~ai1 (|2006f ). 



sitive to the hardness of the photoionizing spec- 
trum. Comparing to the Cloudy model grid of 
Voit (1992), the neon line ratios are consistent 
with values of the ionization parameter, log t/, in 
the range —2 to —2.5 and, assuming a power-law 
photoionizing spectrum oc v a , power-law indices, 
a, in the range —1 to —1.5. 

FigurelUplots the [Nev]A14.3/imvs. [Nem]A15.5 /<m 
line flux densities. The ratio of these lines varies 
stron gly with ionization parameter (ex U 2 , IVoit I 
1992). The strong correlation between these line 
flux densities shown in Figure |4] demonstrates that 
there is little variation in ionization parameter 
amongst the 3C objects in which they are both 
detected. 

6.2. Star-formation Contribution 

The observed infrared emission from AGN 
can have contributions from thermal and non- 
thermal processes. As well as emission from 
warm dust surrounding the nucleus, which peaks 
in the mid-infrared, dust at lower temperatures 
heated by star-formation ca n produce emission 
that peaks in the far- infrared. ISturm et al. ( 2002) 
demonstrated the use of the [Ne v]A14.3/tm to 
[Ne n] A12.8 /tm ratio to determine the AGN contri- 
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Fig. 6.— The thermal fraction, /^°™ = 
Msum / Mlum f° r galaxies (unfilled histogram) and 
quasars (cross-hatched histogram). The BLRG 
3C 22 is included with the quasars. 



bution in composite AGN/starburst objects. For 
our sample, however, we can only constrain the 
[Nev]A14.3/zm to [Ne ill A12. 8 <jm ratio in a smal l 
number of cases. Recently, iBrandl et all (|2006h 
have shown that the 15/30 /im flux density ratio 
is also a very good AGN/starburst discriminator. 
As shown in Figure O almost all objects in the pri- 
mary sample exceed the threshold 15/30 /im flux 
density ratio (~ 0.2) at which the AGN contribu- 
tion dominates (jBrandl et alj |2006). We therefore 
conclude that star-formation does not contribute 
significantly to the measured infrared emission for 
the majority of objects in the primary sample. 



6.3. Non-thermal Contribution 

Synchrotron emission from optically thin radio 
lobes has flux-density spectral indices a « —0.7, 
so the contribution to the infrared is typically 
well over an order of magnitude below the ob- 
served flux density. However, Doppler-boosted 
synchrotron emission from dense, compact regions 
is a potentially serious contaminant of the infrared 
emission. The relative contribution of thermal 
and non-thermal processes to the infrared emis- 
sion of AGN has been th e subject of many stud- 
Sanders et al. ( 19891 ) were among the first to 
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advance thermal dust emission as the dominant 
infrared emission process in quasars, based on a 
study of Palomar-Green (PG) quasars. However, 
claims that quasars are brighter sources of infrared 



emiss i on th a n radio ga l axies (e.g.. lHeckman et al 
1991 11994k iHes et al.l Il995f ) have provided evi- 
dence for a beamed component in the infrared 
emission of FRII sources. Based on a spec- 
tral decomposition for three 3C quasars (3C47, 
3C 207 and 3C 334 ) with IRAS 60 /mi detections, 
van Bemmel et all (1998) found evidence for a 
non-thermal contribution to the IRAS flux den- 
sity in one object (3 C 207). An orien t ation- based 
model described by iHoekstra et al. (1997) was 
consistent with a significant non-thermal contribu- 
tion to the IRAS 60 /um emission of FRII sources 
with higher core dominance. 

The Spitzer data presented here have provided 
further constraints on the role of non-thermal pro- 
cesses in the infrared emission of radio galaxies and 
quasars. In § 15. 2\ we fitted the SEDs of the ob- 
jects in our sample with broadband spectral com- 
ponents representing synchrotron emission from 
the radio lobes and jet as well as thermal emis- 
sion from circumnuclear dust (see Fig. [3]). In this 
way, the non-thermal contribution to the 15 and 
30 /im luminosity was estimated, and the correc- 
tions applied to the data are given in Table[5]as the 
thermal fraction, / thcrm = Lthemy^total at 15 and 
30 fim. Figure [5] shows a histogram of the thermal 
fractions at 15 fim for the quasars and galaxies in 
the primary sample. For the galaxies, the thermal 
contribution is always > 80%; for the quasars, the 
thermal contribution is in the range 10-100%. All 
objects, quasars and galaxies, were fitted in a con- 
sistent manner; however it is only the quasars in 
the sample which were estimated to have a non- 
thermal contribution > 20%. 

Since the errors on the estimated thermal frac- 
tions are likely dominated by model-dependent un- 
certainties, as a consistency check we also u s ed the 
orientation-based model of IHoekstra et al" (1997) 



to predict the non-thermal contrib ution as a func- 



tion o f core dominance. Following Hoekstra et al 



(1997), we fit the ratio of rest-frame 30 /im to 178 
MHz flux density for the primary sample with an 
equation of the form, 



ptotal 
r 30 fan 



78 MHz 



— B DsQ/im (R + Cso/im), (5) 



Fi 

where, F 30tinl and Fn 8 MHz are the rest-frame 
30 /im and 178 MHz flux densities, B is the 
ratio of the rest-frame extended emission at 
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Fig. 7. — Luminosity (WHz _1 sr _1 ) at rest wavelengths 15/zm (left) and 30 pm (right) vs. luminosity at 
rest-frame 178 MHz for objects in the primary sample, before (top) and after (bottom) the correction for 
non-thermal emission has been applied. Symbols are as in Figure [TJ Objects with high thermal infrared 
luminosity for their radio-lobe power are labeled. 
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Table 6 

i?dr STATISTICS FOR QUASARS AND GALAXIES. 







Quasars 


Galaxies 


Q/G a 


(log R dr (15/im)>.... 
(log fi dr (15/W)).... 
(log fi dr (15pm)).... 
(logR dr (15/im)).... 


total 
thermal 
screen 
mixed 


1.78 ± 0.07 
1.44 ± 0.10 
1.66 ± 0.14 
2.00 ± 0.25 


1.21 ± 0.06 
1.19 ± 0.05 
1.66 ± 0.11 
2.74 ± 0.36 


y + 0.9 
°- 1 -0.7 

1 o + ' 5 
1 - 8 -0.4 

1 +0 ' 5 

^'"-0.3 

2 +0 3 


Mcd(log _R dl (15Mm)) 
Mcd(log _R d r(15Aim)) 
Mcd(log _R d ,-(15A<m)) 
Mcd(log R dl (15tim)) 


total 
thermal 
screen 
mixed 


1.85 ± 0.07 
1.42 ± 0.10 
1.67 ± 0.14 
1.72 ± 0.26 


1.29 ± 0.06 
1.21 ± 0.05 
1.58 ± 0.11 
2.41 ± 0.37 


3 K + 0-9 
•>-°-0.7 

1 6 +0 5 

i n + 0.6 
1 - z -0.4 

2 +0 ' 4 


(log_R dr (30Aim)).... 
(log_R. dr (30Ami)).... 
(log.Rdr(30Ami)) .... 
(log B dT (3Q^m)) .... 


total 
thermal 
screen 
mixed 


1.63 ± 0.08 
1.33 ± 0.09 
1.49 ± 0.11 
1.74 ± 0.16 


1.06 ± 0.04 b 
1.06 ± 0.04 b 
1.15 ± 0.10 b 
1.52 ± 0.22 b 


y+0.8 
' 3 - ' -0.7 

1 9 + °' 5 
9 9 + O.9 
z - -0.6 

1 7+ 14 

1 - 1 -0.8 


Mcd(logi?. d r(30/im)) 
Med(Iogiidr(30/im)) 
Med(Iogiidr(30Atm)) 
Med(logiidr(30Atm)) 


total 
thermal 
screen 
mixed 


1.62 ± 0.08 
1.26 ± 0.10 
1.37 ± 0.11 
1.50 ± 0.18 


1.19 ± 0.04 c 

1.20 ± 0.04 c 
1.50 ± 0.06 c 
1.95 ± 0.22 c 


9 y+0.6 
z - ' -0.5 
i 9 + 0.3 
^ -0.3 

7 +0 ' 2 
4+ ' 3 



Note. — The table gives the mean and median of log Rt± r — 
\og(vL u (IR) / vX^ u (178 MHz)) , at the characteristic rest wavelengths of 15 and 
30 /i-m, for quasars and galaxies in the sample with estimates of the infrared 
luminosity and dust optical depth. Also tabulated is the ratio of the mean 
or median R^ T for quasars to that for galaxies. In each case, the statistic is 
computed for the 'total' Rd T (i.e. as derived from the data), the 'thermal' 
Rd T (i-e. corrected for non-thermal contributions), and finally corrected for 
both non-thermal contributions and dust extinction, using the 'screen' and 
'mixed' optical depths. The distribution of R^ r is adequately described as 
log-normal. We therefore estimate the mean of R^ r as 10^ logJ?dr \ and the 
resulting standard error is asymmetric. Where the data include upper limits, 
the statistics arc computed using the Kaplan -Meier (KM) estimate from the 
ASURV statistical package ( Isobc ct al. 1986). The 15 fim lower luminosity 
limit for 3C 200 is treated as a detection. The BLRG 3C 22 is included with 
the quasars. The number of quasars and galaxies included in the analysis 
(number of detections in parentheses) at 15 /xm is 16 (16) and 13 (13); at 
30 it is 16 (16) and 13 (2). 

a This column gives the ratio of the estimated mean and median i?dr for 
the quasars to that of the galaxies. 

b This estimate is biased since a censored datum was changed to a detec- 
tion for the KM computation. Since the KM method does not provide an 
estimate of the error on the median, we estimate the error by treating all 
measurements as detections. 

c The KM estimate of the median falls outside the range of observed values 
due to the heavy censoring of data. Instead, wc estimate the median and 
standard error by treating the upper limits as detections 
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Fig. 8. — Luminosity (WHz _1 sr _1 ) at rest wavelengths 15 /xm (left) and 30 /xm (right) vs. luminosity at 
rest-frame 178 MHz for objects in the primary and supplementary samples. Symbols are as in Figure[TJ FRI 
galaxies from the supplementary sample are plotted as open stars. Data for objects from the supplementary 
sample are plotted using smaller symbols. The infrared luminosity of objects in the primary sample has been 
corrected for non-thermal emission, while that from objects in the supplementary sample has not. 



5 GHz and 178 MHz, £>3o Mm is the ratio of the 
beamed emission at 30 /im and 5 GHz, Cso^m = 
Ho^m/ i F i GHz £>30Mm) is the ratio of the isotropic 
30 /mi emission to the product of the extended 
emission at 5 GHz and D 30 ^ m , and R is the core 
dominance parameter. In this model, D is a mea- 
sure of the strength of the beamed component of 
emission and C is a measure of the value of R at 
which the core emission starts to dominate over 
the extended emission. In order to help constrain 
the high core-dominance region, we also included 
the blazars 3C 454.3 from the supplementary sam- 
ple and 3C345 based on IRAS detections in the 
literature. Despite the many upper limits, we per- 
form this fit at the characteristic wavelength of 
30 /xm, rather than 15 /xm, in order to minimise 
the effect of any correlation of optical depth with 
core dominance. Fitting Equation [5] to the data 
(Fig. [HI left), we obtain a value for C^o^m of 0.05, 
from which the non-thermal fraction as a function 
of core dominance can be estimated as, 



rnon- therm 
/ 30 fim 



R 



R + C30/J 



(6) 



Figure [5] (right) shows the non-thermal frac- 
tions estimated f rom the spectral fittin g and the 
predictions of the LHoekstra et al.l ( 19971 ) model are 
plotted as the solid curve. The dashed curve shows 
the predictions based on the 30 /xm flux density 
corrected for the estimated screen optical depth 
(C3o M m = 0.09). The model predicts a signifi- 
cant non-thermal contribution for the objects in 
the primary sample with higher core dominance. 
Although there is considerable scatter, reasonable 
agreement can be seen between the model predic- 
tions and the results of the spectral fitting. The 
broad- lined radio galaxy, 3C 22, has a significantly 
higher non-thermal estimate from the spectral fit- 
ting than predicted by the model. However, this 
object also has a F 30tim /Fn 8 MHz ratio more sim- 
ilar to the quasars at higher core dominance. The 
quasar 3C 380 has an estimated non-thermal frac- 
tion of around 0.1 from the spectral fits although 
it has the highest core dominance of the primary 
sample. However, this object is an outlier with a 
low i ; 3o At m/-Fi78MHz ratio. We note that the fitted 
value for C3o pm of 0.05-0.0 9 is much less than tha t 
for Ceo^m of 2.6 found by iHoekstra et all ( 1997 ). 
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indicating that the onset of core-dominated emis- 
sion at 30 (im in the primary sample occurs at 
larger inclinations from the radio axis than found 



rgei_ 

by iHoekstra et al.l (|1997T ) at 60 /jm for their sam- 
ple of IRAS-detected 3C and 4C objects. 

Figure [7] shows the 15 and 30 /im luminosity vs. 
178 MHz luminosity for objects in our sample, be- 
fore and after the correction for non-thermal emis- 
sion is applied. Figure [5] shows the non-thermal 
corrected luminosities, but also includes objects 
from the supplementary sample. We were un- 
able to fit dust models for 3C 427.1 (labeled in 
Figs. [7] and HJ since no IRS or MIPS detections 
were obtained for this object. However, we have 
assumed a thermal fraction of 1.0 since it is a lobe- 
dominated galaxy. Comparing the 15 /im luminos- 
ity of the objects in the primary sample before and 
after (Fig. [7]) the non-thermal correction has been 
applied, we can see that the effect of the correc- 
tion is to reduce the luminosity of most quasars to 
a level comparable with that of the galaxies. Note 
that any apparent correlations in the luminosity- 
luminosity plots presented here may be due to a 
mutual dependence on a third variable (e.g., red- 
shift). 

Considering the primary and supplementary 
samples together, we also note that there appears 
to be an upper envelope to the 15 /im luminosity 
of w 10 24 - 5 \VHz~ 1 sr" 1 (see Fig. E). Of the ob- 
jects in the primary and supplementary samples, 
only the blazar, 3C 454.3 has a luminosity greater 
than 10 25 WHz -1 sr _1 . We will discuss this upper 
envelope further in the next section. 

6.4. Comparative Luminosity of 3CRR 
Radio Galaxies and Quasars 

In order to compare the infrared luminosities of 
quasars and radio galaxies we must first account 
for the variation of central engine power amongst 
objects in our sample. To do this, we normalise the 
infrared emission by the i sotropic low-frequenc y 
radio luminosity (see e.g., iHeckman et all 1 19921 ). 
Figure [TD] shows i? dr = ^L„(IR)/i/L iy (178MHz) 
for our sample at 15 and 30 /im. Figure fTTI is the 
same, but includes objects from the supplemen- 
tary sample. Table 6 shows the mean and median 
logi?dr for quasars and galaxies in the primary 
sample, as well as the ratio between the mean or 
median i?dr for quasars to that for galaxies. For 
the purposes of this comparison, the BLRG 3C 22 



has been included with the quasars. Only objects 
with characteristic luminosity and optical depth 
estimates are included in the calculation of these 
statistics. Due to the low detection rate for galax- 
ies at 70 /im using MIPS, we base our conclusions 
in this section on the 15 /im statistics. 

The mean i?dr for quasars is ~ 4 times higher 
than for galaxies at 15 /im. Subtracting the non- 
thermal emission brings this down to k 2. The 
fitted optical depths allow us to go further and 
ask whether obscuration by dust can account for 
this residual anisotropy of the thermal emission 
between quasars and galaxies. Using the fit- 
ted screen optical depths to correct for dust ex- 
tinction, we find that the thermal emission from 
quasars and galaxies is on average equal at 15 /im. 
The mixed optical depths would imply intrinsi- 
cally brighter galaxies than quasars. Since there is 
no evidence that galaxies are brighter than quasars 
in the far infrared, where the torus is expected to 
be optically thin, we conclude that the actual dis- 
tribution of warm and cool dust is better approx- 
imated by the screen model. Therefore, around 
half of the difference in mid-infrared brightness 
between quasars and galaxies can be attributed to 
Doppler-boosted synchrotron emission present in 
some quasars, but not in galaxies. The other half 
is consistent with dust absorption in the galaxies 
but not in the quasars. 



Haas et al.l (|2004l ) find no difference in the dis- 



tribution of i?dr for quasars and galaxies for their 
sample of 3CR radio galaxies and quasars; how- 
ever, we believe that this is not inconsistent with 
the re sults presented he re for the following rea- 



Haas et al.l (|2004|) exclude from their sub- 



sample those objects which do not exhibit a promi- 
nent thermal bump, thereby excluding objects 
with significant non-thermal contributions. Also, 
their mean i?dr is computed at rest-frame 70 /im 
where the effect of optical depth is much less than 
at 15 /im. 

Figure [12] plots the screen, Tg c j^ m , and mixed, 
Tg "* m , optical depths against the core dominance 
parameter, R = F COTO / intended (see Table 7). Our 
dust model (see § 15. 2|) does not, a priori, impose 
any asymmetry in the distribution of dust sur- 
rounding the nucleus. From Figure [12] there is evi- 
dence of an anti-correlation between optical depth 
and core dominance, from which we can infer an 
equatorial distribution of dust. The median screen 
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Fig. 9. — The orientation-based model (solid curve) of iHoekstra et al.l (|1997l ) is fitted to the ratio of the 
rest-frame 30 /im to 178 MHz flux density (left). Based on the fit to the data, the model can be used to 
predict the non-thermal fraction as a function of orientation. We plot the non-thermal fraction for each 
source in the primary sample estimated from the spectral fitting (right) and also the model prediction (solid 
curve). The dashed curve shows the predictions based on the 30 /jm flux density corrected for the estimated 
screen optical depth. The model predicts significant non-thermal contributions for the more core dominated 
objects in the primary sample. Although there is considerable scatter, reasonable agreement can be seen 
between the model predictions and the results of the spectral fitting. Symbols are as in Figure [TJ 
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Fig. 11.— i? dr (15/mi) (left) and i? d r(30/im) (right) vs. rest-frame 178 MHz luminosity (WHz^'sr 1 ) for 
objects in the primary and supplementary samples. Objects in the primary sample have been corrected for 
non-thermal emission. Symbols are as in Figure [TJ Data for objects from the supplementary sample are 
plotted using smaller symbols. An upper envelope, defined by setting i?d r (15^m) = K x (10 24 - 5 /^i78MHz), 
where K = 3 x 10 8 /(15 x 178), is indicated by the dashed line. The blazar, 3C 454.3, lies above this upper 
envelope. 



(mixed) optical depth for objects with R > 10 2 
is 0.4 (1.0) and with R < 10~ 2 is 1.1 (3.0). This 
is in line with the predictions of orientation-based 
unification of 3CRR quasars and radio galaxies. 

At lower radio luminosities than those of the 
primary sample, we can see from Figure [ITJ that 
there are quasars and galaxies of similar radio lu- 
minosity but separated by over three orders of 
magnitude in R^ T , suggesting that differences may 
exist between these objects which cann ot be ex 



plain e d by orientation alone (see e.g., IShi et al 



20051 : lOgle et all [2006). The lack of similar ob- 



jects at the greater radio luminosities spanned by 
the primary sample may be due to selection ef- 
fects, since the most radio-luminous objects in the 
flux-limited 3C sample are farther away. There 
is, however, one notable exception: the galaxy 
3C 427.1, which was not detected using MIPS or 
IRS. Since the MIPS 24 /im measurement for this 
object is at ss 15 /im in the rest-frame, we ob- 
tain logi?dr(15/^m) < (see Fig. [Til left), which 
is over two orders of magnitude below quasars of 
similar radio luminosity. This galaxy may there- 
fore be different to the quasars and other galaxies 
in the primary sample, perhaps with an intrinsi- 
cally weak source of optical/UV compared to other 
powerful radio sources. 



We note that there is an apparent upper enve- 
lope in i?dr(15 /im), due to the previously noted 
upper envelope of Li 5 ^ m < 10 24 5 WHz _1 sr _1 . 
This is shown by the dashed line in Figure [TJJ 
(left), which is defined as i?dr(15/im) = K x 
(10 24 - 5 /ii 78 MHz), where K = 3 x 10 8 /(15 x 178). 
The blazar 3C 45 4.3 lies above this envelope. 
lHaas et al] (2004) point out that a decline of 
i?dr with radio luminosity is expected from the 



"receding-torus" model of lLawrence I (|1991l ). yet 



find no evidence for such a decline in their sam- 
ple of 3CR objects and so suggest a modification 
to the original proposal in which the torus scale 
height remains constant with luminosity. Taking 
all quasars in the primary and supplementary sam- 
ples (except the blazar 3C 454.3) and computing 
the mean (non-thermal corrected) logi?dr(15/zm) 
in two bins defined by their median radio lumi- 
nosity (« 10 27 - 5 WHz^sr 1 ), we find that R Al 
for quasars is around 4^ times greater in the 
lower radio luminosity bin. For galaxies, i?dr is 
approximately equal in both bins. The observed 
decline in i?dr for quasars is of low significance; 
however if the observed upper envelope in 15 /im 
luminosity is shown to apply at higher radio-lobe 
luminosities than those of the Spitzer sample, then 
this would provide further support for the receding 



21 



torus model as originally formulated bv lLawrence 
(jl99lh . 



7. CONCLUSIONS 

We have observed a sample of extremely pow- 
erful radio sources with the IRS and MIPS instru- 
ments aboard Spitzer. The sample was selected 
at low frequency, and is orientation independent. 
Our aim was (i) to determine the mid-infrared 
properties of an orientation- unbiased sample of the 
most powerful radio sources, (ii) to assess the rel- 
ative contribution of the AGN, star-formation and 
Doppler-boosted jet synchrotron and (iii) to study 
the relative importance of physical differences and 
orientation effects in powerful radio sources. 

Our main results can be summarised as follows. 

- We have detected extremely luminous emis- 
sion in the mid-infrared for our sample of ra- 
dio galaxies and quasars. All but one object 
observed in our sample with MIPS were de- 
tected at 24 /im, with specific luminosities, L u > 
lO^WHz^sr" 1 . 

- The MIPS and, in particular, IRS measure- 
ments provide powerful constraints on the spectral 
energy distributions of the sources in our sample. 
They have allowed us to fit continuum synchrotron 
and dust components representing the dominant 
emission processes and thereby estimate the non- 
thermal contribution to the mid-infrared luminos- 
ity. We find that non-thermal processes can con- 
tribute a significant proportion (up to 90%) of the 
measured infrared emission in some quasars. 

- Warm, optically thin dust components with 
temperatures in the approximate range 100-300 
K, obscured by cooler, optically thick dust, pro- 
vide very good fits to the IRS spectra (when com- 
bined with appropriate synchrotron components 
from the lobes and jet). The optical depths re- 
sulting from the fits are typically tJ" < 2 
(Ay < 40) when using a dust screen model and 
r 97^m ~ 8 (^v ^ 150) when considering a mixed 
extinction scenario. 

- For the most powerful radio sources (Lng, mhz ~ 
10 26 5 WHz -1 sr _1 ), the mid-infrared luminosity 
appears to have an upper envelope of around 
Li5p, m = 10 24 - 5 WHz -1 sr _1 . If this envelope ap- 
plies at even higher radio-lobe luminosities than 
those of the Spitzer sample, this may p rovide sup- 
port for the "receding torus" model of lLawrence 



1)199111 . 

- Normalising by the isotropic radio luminosity 
(to account for variations in central engine power), 
the quasars are ~ 4 times brighter than the galax- 
ies in our sample at 15 /Ltm. Half of this difference 
can be attributed to the presence of non-thermal 
emission in the quasars but not in the galaxies. 
The other half is consistent with absorption by 
a dust screen in the radio galaxies but not the 
quasars. 

- We find an anti-correlation between the es- 
timated optical depth and core dominance. The 
median screen optical depth for objects with core- 
dominance factor R > 10~ 2 is « 0.4; for objects 
with R < 10~ 2 , it is ps 1.1. This is consistent 
with the picture of the MIR-emitting dust in FR II 
galaxies being observed through an equatorial dis- 
tribution of cooler dust, while for the quasars, the 
MIR-emitting dust is seen directly as the torus is 
viewed face-on. 

Previous investigators have shown that quasars 
are brighter than galaxies in the mid-infrared, in 
agreement with the results presented above. It is 
also clear from previous work that beamed syn- 
chrotron emission and dust extinction modulate 
the mid-infrared emission of quasars and radio 
galaxies to some degree. In the work presented 
here, however, we have quantified these effects for 
an orientation-unbiased sample of powerful radio 
sources and shown that once they are taken into 
account, quasars and radio galaxies are on average 
equally bright in the mid-infrared. 
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Laboratory, California Institute of Technology un- 
der a contract with NASA. This research has made 
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Laboratory, California Institute of Technology, 
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Fig. 12. — The 9.7/zm screen (left) and mixed (right) optical depths estimated from from the fits to the data 
versus the core dominance parameter, R — F coro /F ox t C ndcd- The right-hand ?/-a xis on each plo t indicates the 
visual extinction, based on Ay/Tg^m — 18.5, typical for the local diffuse ISM (jDraine 1 12003). Symbols are 
as in FigurefT] The error bars indicate the formal uncertainties in the fits. From the apparent anti-correlation 
of optical depth with core dominance, we infer an equatorial distribution of cool dust consistent with the 
"dusty torus" of orientation-based unification schemes. 
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B. FITS TO THE SPECTRAL ENERGY DISTRIBUTIONS 




Fig. 13. — Fits to the IRS spectra and photometric data (left). The fit to the IRS data (where available) is 
also shown in detail (right). Filled circles and triangles represent detections and upper limits respectively. 
Symbols within circles represent the MIPS measurements. For clarity, the MIPS 24 /im measurement is 
excluded from the plot unless no IRS spectrum was obtained. Data from the literature which appear 
discrepant with Spitzer data are plotted as open symbols and excluded from the fits. The characteristic 
30 /im flux density is plotted as an open circle within a square. The dark-shaded error contours represent the 
uncertainty in the IRS spectrum obtained by differencing the two nod positions. The light-shaded contours 
represent the formal la uncertainty in the fitted models (note that the fit uncertainties are greatest in regions 
where the dust and synchrotron components overlap). IRS data with signal-to-noise ratio < 1 are connected 
with dotted lines and excluded from the fit. 
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Fig. 14. — The fit to the IRS spectrum for three quasars with warm dust + lobe synchrotron + cool dust 
(left) and warm dust + lobe synchrotron + jet synchrotron (right). The importance of the constraints from 
the IRS spectrum on the fits is clear. 
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C. NOTES ON THE FITS FOR SOME INDIVIDUAL OBJECTS 



3C6.1: IRS SL2 was not detected for this source and SL1 alone provides poor constraints on the optical 
depth from the fits. 

3C 22: The best fitting combination of components includes a jet functional form, since this component 
greatly improves the fit to the IRS spectrum. The next-best fit involved two dust components and a lobe 
synchrotron component, as shown in Figure [T4"l but the fit to the IRS spectrum is much poorer. Since this is 
a broad-line galaxy, it is not unreasonable to expect some contribution from the jet, however sub-millimeter 
measurements of this source are required to settle the issue. 

3C 200: No IRS spectrum was available for this galaxy due to an IRS peak-up failure. Instead, a warm dust 
component with the temperature fixed to T = 185 K was fitted to the MIPS 24 /im detection. 
3C 207: This object has been examined in detail bv Ivan Bemmel et al. (1998) who performed a spectral 
decomposition based on simultaneous cm and mm data and estimate that substantia l non-thermal c o ntam - 
ination of the IRAS 60 /im flux density is possible. The orientation-based model of iHoekstra et al.l (jl997l ) 
yields an estimate for the non-thermal contribution of around 15% at 60 /im. We estimate a non-thermal 
contribution of approximately 90% at rest-frame 30 /im. 

3C 216: This is a quasar with known super luminal motion and beamed jet emission is clearly present in the 
radio and millimetric data. A combination of three components best fits the data: a dust component with 
T = 158 K, low-frequency synchrotron emission from the lobes, and Doppler-boosted synchrotron emission 
from the jet. We estimate that the jet contributes around 90% of the flux density at 15 /im. 
3C 220.3: This galaxy provided the only MIPS detection at 160 /im, as well as detections at 24 and 70 /im. 
The measured 160 /im flux density was very high, at 136 mJy No IRS spectrum was available for this galaxy 
due to an IRS peak-up failure. Instead, a warm dust component with fixed temperature T = 185 K and a 
cool dust component whose temperature was a free parameter were included in the fits. The resulting best- 
fitting cool component has a temperature T = 31 K and comprises 69% of the total fit ted dust luminosity, 
a proportion comparable to that seen in starburst galaxies (see e.g., Marshall et al. 2006f) . We omitted a mm 
upper limit from the fit in order to achieve a good fit to the radio data using the lobe synchrotron functional 
form. 

3C 272: No IRS spectrum was available for this galaxy due to an IRS peak-up failure. Instead, a warm dust 
component with the temperature fixed to T = 185 K was fitted to the MIPS 24 and 70 /im detections. 
3C 280: We have not included an apparently discrepant ISO detection from I Andreani et al.l (|20021 ) in the 
fit. 

3 C 286: We have not included apparently discrepant ISO detections from iMeisenheimer et al.l (|2001l ) and 



Andreani et all (|2002l ) in the fit. 



3C 289: No IRS spectrum was available for this galaxy due to an IRS peak-up failure. Instead, a warm dust 

component with the temperature fixed to T = 185 K was fitted to the MIPS 24 and 70 /im detections. 

3C 309.1: The best fit for this SSC quasar has synchrotron components from the radio lobes and jet, as well 



as a d ust component with T = 181 K. However, there is a discrepant ISO detection at 60 /im (jAndreani et al 



2002) which we have excluded from the fi t. We were unab le to achieve a good fit to the IRS data without 



also excluding the upper limit at 850 /im (|Haas et al.l 120041 ) . 



3C 325: The ISO 60 and 90 /im detections from lAndreani et alJ (|2002l ) appear inconsistent with MIPS upper 
limits and so they have been excluded from the fit. 

3C 334: As well as lobe synchrotron and warm dust, the best fit includes a jet functional form which exceeds 
a sub-millimeter upper limit. Although some far- infrared data are not well fit by the model, a cool dust 
component was not also included in the best fit since the addition of this component resulted in a confidence 
level for the F-test statistic of just 40%. The inclusion of a jet in the best fit is largely due to the IRS 
data (see Fig. [T4|) and we estimate the non-thermal c ontribution to be approx imately 20% at rest-frame 
30 /im. This object has also been examined in detail by van Bemmel et all (|l998l ) who performed a spectral 
decomposition based on simultaneous cm and mm data and estimate a negligible non-thermal contribution 
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to the IRAS 60 /xm flux density. In addition, the orientation-based model of iHoekstra et al. (Il997 ) yields a 
non-thermal contribution of much less than 15% at 60 /jm. However, as noted bv lvan Bemmel et al 
the existence of a variable sub-millimeter component, not accounted for in their analysis, is likely in this 
superluminal object. 

3C 427.1: No constraints were available in the infrared since this galaxy was not detected with MIPS or 
IRS. 
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1.0 


3C272 


G 


6.9 





14 


1.0 


3.0 










185 










6.6 


1.0 


1.0 


3C280 


G 


2.3 





12 


1.0 


3.0 


0.0 


300.0 


252.1 


80.1 


310 


1.5 ± 


1 


7.9 ± 


4 


31.6 


0.9 


0.8 


3C289 


G 


3.3 





14 


1.0 


3.0 










185 










12.2 


1.0 


1.0 


3C292 


G 


0.8 





16 


6.2 


3.0 










283 


1.3 ± 


4 


3.5 ± 1 


3 


2.3 


1.0 


1.0 


3C325 


G 


2.2 





18 


7.4 


3.0 


0.0 


3.0 


193.5 




172 


0.0 ± 


1 


0.1 ± 


2 


14.4 


0.8 


0.9 


3C330 


G 


6.7 





12 


1.0 


0.1 








80.1 


293 


1.2 ± 


1 


4.0 ± 


4 


5.9 


1.0 


1.0 


3C337 


G 


3.7 





13 


1.0 


0.3 










310 


1.1 ± 


3 


2.9 ± 2 


6 


1.1 


1.0 


1.0 


3C340 


G 


1.5 





15 


2.7 


3.0 










275 


1.6 ± 


2 


5.7 ± 1 


3 


3.4 


1.0 


1.0 


3C 343.1 


G(SSC) 


5.5 





62 


307.4 


3.0 


-1.2 


9.4 


600.0 




190 


0.4 ± 


3 


1.0 ± 


8 


2.1 


0.9 


0.9 


3C352 


G 


2.6 





16 


1.0 


3.0 


-1.0 


20.7 


60.0 




222 


0.7 ± 


3 


1.7 ± 


9 


2.0 


1.0 


0.9 


3C 427.1 


G 


9.6 





15 


1.0 


3.0 


























3C441 


G 


3.2 





13 


1.0 


3.0 










292 


2.0 ± 


3 


6.6 ± 2 


2 


2.7 


1.0 


1.0 



Note. — The table gives the free parameters for the synchrotron fits, /3, v± , ^i, obc , c*2 > y 0i and defined in the lobe synchrotron fitting function, log F u (is) oc — {3 (log is — 

log ft ) 2 + log[e — u c ) ] and the jet synchrotron fitting function, F u (i/) a e~ T ^' (fo/^} -01 (\. — e~^ u o/ 1 ') 1 2 ^ & ~( v i v c ) _ Also tabulated arc the reduced chi-squarcd, i 

the temperatures of the best-fitting cold and warm components, T co ^ and T warm , the screen and mixed optical depths, r scr and T mix , the integrated luminosity of the fitted dust 
components, L(j ust and the ratio of thermal to total luminosity at 15 /jm and 30 /Am, fj^ 1 ^ 1 and f^rwi™ ■ 
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Table 7. Characteristic luminosities for the primary 3CRR Sample. 



Name 


Type 


^178/(l+z) MHz 

[Jy] 


F 5/(l + z) GHz 

[Jy] 


F 5GHz 
[mjy] 


F 15(l + z)^m 

[mJy] 


F 30(l + z)nm 

[mjy] 


log L 17s M Hz 
[WHz -1 sr" 1 ] 


log£l5nm 
[WHz -1 sr _1 ] 


[WHz -1 sr 


3C 138 


Q(SSC)(SL) 


28.9 


6.0 


94.0 


18.9 


21.7 


27.5 


24.4 


24.4 


3C 147 


Q(SSC)(SL) 


74.5 


11.8 


2500.0 


20.1 


39.9 


27.7 


24.1 


24.4 


3C 175 


Q 


31.0 


1.0 


23.5 


18.2 


22.2 


27.6 


24.3 


24.4 


3C 196 


Q 


118.0 


7.3 


7.0 


22.7 


13.2 


28.3 


24.6 


24.3 


3C 207 


Q 


111 


1.8 


510.0 


12.1 


17.7 


27.0 


24.1 


24.2 


3C 216 


Q(SL) 


34.9 


1.9 


1050.0 


22.0 


47.7 


27.5 


24.3 


24.6 


3C 254 


Q 


41.4 


1.0 


19.0 


13.8 


11.9 


27.7 


24.2 


24.1 


3C 263 


Q(SL) 


21.3 


1.6 


157.0 


27.7 


29.7 


27.3 


24.4 


24.4 


3C 275.1 


Q 


26.3 


1.1 


130.0 


8.1 


15.8 


27.2 


23.7 


24.0 


3C 286 


Q 


38.3 


10.6 




9.7 


21.9 


27.8 


24.2 


24.5 


3C 309.1 


Q(SSC) 


33.2 


5.4 


2350.0 


22.8 


31.6 


27.8 


24.6 


24.7 


3C 334 


Q(SL) 


14.0 


0.8 


111.0 


33.0 


51.6 


26.9 


24.3 


24.5 


3C 336 


Q(SL) 


19.9 


1.5 


20.4 


5.3 


10.1 


27.6 


24.0 


24.3 


3C 343 


Q(SSC) 


19.2 


3.0 


<300.0 


16.6 


32.1 


27.6 


24.5 


24.8 


3C 380 


Q(SSC)(SL) 


92.3 


9.0 


7447.0 


30.4 


32.5 


28.0 


24.5 


24.5 


3C 6.1 


G 


19.7 


2.0 


4.4 


3.8 


<9.7 


27.5 


23.7 


<24.2 


3C 22 


G 


20.3 


1.4 


7.8 


16.7 


20.8 


27.6 


24.5 


24.6 


3C 184 


G 


22.0 


1.4 


<0.2 


4.2 


8.7 


27.7 


23.9 


24.3 


3C200 


G 


13.9 


0.9 


35.1 


>1.1 


<4.3 


26.8 


>22.7 


<23.3 


3C 220.1 


G 


20.3 


1.0 


25.0 


2.7 


<5.0 


27.2 


23.3 


<23.6 


3C 220.3 


G 


29.0 


1.2 


<0.2 


5.2 


13.3 


27.5 


23.7 


24.1 


3C 268.1 


G 


32.3 


4.5 


2.0 


3.8 


<5.8 


27.8 


23.9 


<24.1 


3C 272 


G 


12.8 


0.8 




2.8 


7.5 


27.4 


23.7 


24.1 


3C280 


G 


36.4 


3.1 


1.0 


13.1 


16.5 


27.9 


24.4 


24.5 


3C289 


G 


21.6 


1.3 


10. (i 


5.0 


13.6 


27.6 


24.0 


24.4 


3C292 


G 


11.7 


1.1 


1.0 


2.2 


<5.6 


27.1 


23.4 


<23.8 


3C325 


G 


28.1 


2.2 


2.4 


5.5 


<6.7 


27.9 


24.2 


<24.3 


3C330 


G 


38.6 


3.3 


0.7 


7.1 


<10.8 


27.4 


23.6 


<23.8 


3C337 


G 


20.1 


1.5 


0.3 


1.0 


<3.9 


27.2 


22.9 


<23.5 


3C 340 


G 


15.8 


1.2 


1.1 


2.7 


<6.6 


27.3 


23.5 


<23.9 


3C 343.1 


G(SSC) 


13.7 


1.8 


<200.0 


1.5 


<5.4 


27.2 


23.3 


<23.8 


3C 352 


G 


21.8 


0.9 


3.2 


1.7 


<4.2 


27.5 


23.4 


<23.8 


3C 427.1 


G 


35.0 


1.7 


0.8 


<0.4 




27.4 


<22.4 




3C441 


G 


17.5 


1.4 


3.5 


2.3 


<3.9 


27.3 


23.4 


<23.6 



Note. — The tabic gives the flux density in the observed frame, F^, , at frequencies, v, of 178/(1 + z) MHz and 5/(1 + z) GHz, using the best fitting lobe synchrotron 
model fits (see Table 5) to the radio data. The observed core flux density at a frequency of 5 GHz, F§ < qii z , is also given. Assuming a core flux density spectral index 
a — 0, the core-dominance parameter is calculated as R — F§ *qhz /-^*5/ (1+z) MHz ■ The observed flux density at a wavelength corresponding to 15 /im in the rest-frame, 
^15(1-1-2)^111' * s computed from the IRS spectrum as the median flux density in a 0.2 fim bin centered on 15(l+z) //.m. The flux density at a wavelength corresponding 
to 30 /Am in the rest-frame, ^3o(i-|-z) /im; * s computed from the MIPS 24 and 70 fJ. m flux densities as ^30(1-1-2) = ^24^m(30(l + z )/24) — a , where a is the 24— 70 /Am 
flux density spectral index (F v oc Note that the 30 /im luminosities arc derived from MIPS photometric measurements with a specific spectral response, whereas 

the 15 /tm flux densities arc derived from spectroscopic measurements averaged over a narrow wavelength bin. Care must therefore be taken if directly comparing these 
quantities. The rest-frame luminosities, L\7S,~m.Hz^ -^lS^im and L^^ m , are calculated per steradian as L Uc — ^*"i^ c / ( 1 + s ) -^L / ( ~f~ z )i where v c is the emitted frequency 
and .Dl 1S the luminosity distance. 



